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A perfecdy conducting prolate spheroidal antenna immersed in a conducting medium is analyzed. 
The dimensions of the spheroid are assumed to be small compared to a skin depth in the medium. A 
gain factor for "the antenna as a receiver is computed as a function of the spheroid's axis ratio. The 
prolate spheroid is then assumed to be excited by a circular filament of magnetic current located in 
the plane of the spheroid's minor axes to approximate excitation by a magnetic toroid. Expressions 
for the antenna's input admittance and effective length are derived. Some experimental results veri- 
fying the derivation of the input conductance are included. A computation shows that at VLF, the 
antenna has a ^ of 0.5 or less, and a relatively constant effective length. 



1. Introduction 

It is well known that the receiving properties of a 
small electric loop antenna can be enhanced by loading 
the loop with a permeable core. A theory developed 
by Wait [1953a, 1953b], valid for small loops in a 
uniform field, shows that the permeable core tends to 
concentrate magnetic flux within itself thereby in- 
creasing the effective area of the loop wound around 
the core. An analogous situation occurs with a sphe- 
roidal piece of metal in a medium where conduction 
currents dominate displacement currents. Lines of 
electric current tend to concentrate in the metal which 
is assumed to have a much greater conductivity than 
the dissipative medium surrounding it. W. L. Ander- 
son [1961a, 1961b] seems to have been the first person 
to consider seriously the use of this phenomenon to 
improve the reception of antennas buried in the earth. 
In this paper, the basic ideas of Anderson are modified 
and extended to analyze the prolate spheroidal antenna 
immersed in a conducting medium. 

The general problem of a small electric type antenna 
in a dissipative medium was investigated first by C. T. 
Tai, in which he analyzed the fundamental electric 
dipole in an insulating radome [Tai, 1947], and the thin 
biconical antenna in a spherical insulating radome 
[Tai, 19491. Wheeler [1958] and Hansen [1963] sum- 
marize some properties of the electric type antenna 
and compare these with the properties of the magnetic 
loop antenna in an insulating radome. One purpose 
of the insulating radome in these analyses is to sepa- 
rate idealized circuit theory type sources from the 
dissipative medium. For, if these sources are in con- 
tact with a dissipative medium, then losses become 
unbounded. Equally vaUd methods of ehminating 
these unbounded losses can be made other ways. 



One example is given by King, Harrison, and Denton 
[1961], who consider a small uninsulated cylindrical 
antenna in a dissipative medium. Their method, 
however, requires the current distribution to vanish 
at the ends of the cylindrical antenna, which means 
that their antenna length cannot be too small. Another 
example is given in this paper in which the source is 
specified in terms of a practical toroidal winding. 

Assume a perfectly conducting prolate spheroidal 
antenna that is small with respect to the skin depth 
in the dissipative medium surrounding it. A prolate 
spheroidal coordinate system is used [Magnus and 
Oberhettinger, 1954, p. 149, 157] (see fig. 1). Since 
the spheroid is assumed small, the quasi-static ap- 
proximation to Maxwell's equations is valid in regions 
near the spheroid. Thus, Maxwell's equations are 
taken to be 



VxH = a-E, 
VxE = -M, 



(1) 



where M is the magnetic current density vector. 

In the following sections of this paper, advantage 
is taken of the azimuthal symmetry of the spheriod 
shown in figure 1. With this symmetry it is clear that 
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Figure 1. 



Prolate spheroid in an infinite dissipative medium in 
which cr/(coei) > 1. 
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Ei> = Hv = Hi = and that //»>, Ev, and Ei are independ- 
ent of (f. It then follows that 



VxVxH = 



(tM 



(2) 



reduces to a scalar partial differential equation when 
M is oriented in the (^-direction, and is independent of 
(f. In regions where M = 0, (2) is separable and its 
solutions can be constructed as a product of linear 
combinations of associated Legendre functions^ in 
both variables iq and ^: 



H = J^P]lr^)[AnPUi) + BnQW] ' 



(3) 



2. Receiving Properties 

Let an electromagnetic plane wave be incident on 
the prolate spheroid in figure 1 such that the electric 
vector is parallel to the z-axis. Let the spheroid be 
defined by the surface f = fo. Since the spheroid is 
small, the incident electric field can be assumed uni- 
form in the region about the spheroid. Thus, 

or, in terms of prolate spheroidal coordinates [Magnus 
and Oberhettinger, 1954, p. 149]: 



The second-kind solution in 17, Qliirj), is not used be- 
cause it does not remain bounded when 17 =± 1. 

Using (1), the electric field components are expressed 
as 



and 



Ei=(E')v 



+ E* 



^^-ff Vf2_^2 



PA-n)P\{0, 



Er,-- 



C<J 



V^2-Tj2a^ 



[V^2 _!//], 



=(£')^ 






7,- 



--PA^)n-n)- 



(6) 



Ei = - 



ca 



Ve 



■T7' 



dri 



[Vl-rj^W], 



(4) 



where c is the semifocal length of the prolate spheroid. 
The metrical coefficients for prolate spheroidal co- 
ordinates are listed here for convenience since they 
are referred to often in the remainder of this paper: 



If the metallic spheroid were not present, the current 
enclosed by a circle with radius b located in the 
X, y-plane at the origin would be 



h = 7Th''(jEK 



(7) 



With the spheroid present the current enclosed by the 
same circle would be 



hr) = C 



^2_l 



e-v' 



/l^ = cV(^2-l) (1-7,2). 



(5) 



In terms of the assumptions already stated, the re- 
ceiving properties of the prolate spheroid are derived 
in section 2, the driving point admittance is derived 
and compared with some experimental results in sec- 
tion 3, the effective length and equivalent circuit are 
discussed in section 4, and the conclusions are pre- 
sented in section 5. 



2 Since there can be no confusion, the (^-subscript in H is deleted. Also, no special 
notation for f,/ and Q^ is used to show that |rj| =s 1 and ^ ^ I. In this paper the following 
definitions are used: 



Pl(r)) = 




Pniv) Nl^l, 



piio=ve^^p„i^) 



<?;,(a=VF=lJ^C>„(^) 



^>1. 



For other definitions and properties of associated Legendre functions of the first and second 
kind, see Magnus and Oberhettinger [1954, ch. 4]. 



I = gh. 



(8) 



where ^ is a current collecting gain factor to be derived 
in this section. Anderson [1961b] previously derived 
a similar gain factor using a method different from that 
used here. He assumed a prolate spheroid of finite 
conductivity and then used the formal analogy between 
stationary currents, and static fields in ideal dielectrics 
[Panofsky and Phillips, 1962, p. 120]. For spheroids 
of finite conductivity this method, of course, led to 
electric fields inside the spheroid which are parallel 
to the z-axis. This result is valid only in the strictly 
stationary current case, and does not allow for the 
appreciable skin effect phenomenon in the metal which 
occurs in this quasi-static problem. Although Ander- 
son's gain factor becomes the same as the one derived 
here when the conductivity of his spheroid becomes 
infinite, the method of derivation presented here seems 
preferable. 

The incident field causes induced currents in the 
prolate spheroid and produces a scatter field which 
can be constructed as a weighted sum of Legendre 
functions as shown in (3): 



H^ = ^BnPii7))Qi(^), 



(9) 



where the coefficients An in (3) are set equal to zero 
because Pji(^)—^ '^ as ^^ oo. The coefficients Bn are 
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determined by requiring the total tangential electric 
field at the surface of the prolate spheroid to vanish. 
That is 



^i,(^o,rj) + £-^(|o,r,) = 0. 



(10) 



Using (4) and (9), the only nonzero coefficient is Bi; 
acE'PAM 



5,=-- 



2 <3,(^o)' 



(11) 



The total current flowing through the circle of radius 
b with the conducting prolate spheroid in place is 



/ = 



r2n ri 
Jo Jo 



aE(dr)d(p 



f = fo 



(12) 



where Ei is the total normal field from (4) and (6): 
E> 



Ee=- 



^e-v' 



■P.(r))Pj(^) 






Performing the integrations, one obtains 

PMo) 



I = ttc^(tEW^I-1 
which can be written as 



P\{U)- 



Q<(^o 



Vi-ri'^Pli-n) 

(13) 

Q'Mo) . 

(14) 



where the Wronskian of Pii^o) and ^i(^o) is [Jahnke 
and Emde, 1945, p. 114] 



W[Pr(^o).Qi(^M = 



-1 



(16) 



The radius of the prolate spheroid is b as shown in 
figure 1, and is related to the prolate spheroidal 
surface ^o by 



b=cVa-i- 

Using (7), (14) can now be written as 

7 = ^ . 

iil-Wii^o) 

Therefore the gain defined by (8) is 

1 



(17) 



(18) 



(Cur 



g),R= 20 log. 
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Figure 2. Relative current gain of the prolate spheroidal antenna. 




(a) Prolate spheroid excited by 
a ring of magnetic current. 



(b) Prolate spheroid driven 
by a voltage across an 
infinitesimal gap. 




(a-Wi(^o)- 



(19) 



(c) Prolate spheroid driven 
by a magnetic toroid. 



Fl(;uRE 3. Methods of prolate spheroid excitation. 



3. Input Admittance 

To find the input admittance of a driven prolate 
spheroidal antenna, assume a ring of magnetic cur- 
rent K located at the center of the prolate spheroid to 
be the source, see figure 3a. If the spheroid is sepa- 
rated into halves and driven by a voltage source, then 
the location of the magnetic current ring is on the sur- 
face of the spheroid, (^i = ^o) as shown in figure 3b. 
In case the spheroid is driven by a torroidal winding 
then K is an approximation to the magnetic current of 
the torroid, and is located in the center of the toroid 
cross section, (^i > ^o), as shown in figure 3c. In any 
case 



Figure 2 shows the gain in dB as a function of the axes 
ratio alb of the prolate spheroid. 



''" V Jo" 



Kdl, 



(20) 



767-936 0-65— 8 



1005 



where K is related to the voltage V by 



V=-K, 

and the integration is over the contour specified by the 
ring source. 

In terms of the magnetic current density vector, the 
ring source is defined to be 






(21) 



where dix) represents the Dirac delta-function. The 
primary quasi-static field from the ring source, that is, 
the field produced with no metallic spheroid present, 
can be obtained by substituting (21) into (2), assuming 
a solution of the form specified by (3), applying the 
orthogonality relationship to the functions of 17, and 
applying the method of variation in parameters to the 
resultant differential equation involving f . The result 



IS 



HP = 2H„Pi(v)' 



Q'ni^)Pli^l), 



(?i(^on(a. 



f>6 



(22) 



where 



The field scattered by the perfectly conducting 
prolate spheroid is constructed as follows: 

//^ = 2 CnHnQl(^l)Pi{-nmiO, ^ > f 0, (24) 

^< = l 

where ^0 is the fixed ^ describing the surface of the 
prolate spheroid. The total magnetic field is the sum 
of the primary and scatter fields: 



Ht = HP^H'. 



^^=-iI^^ 



dl. 



(28) 



The expression for Fo in (27) is that part of the total 
admittance Fin which is due to a magnetic toroid alone 
in a conducting medium. The quantity Yq has real 
and imaginary components and has been determined 
by Swain [1965]. The AG expression in (28) is the 
contribution to the total admittance Fjn due to the 
presence of the conducting prolate spheroid. Carry- 
ing out the integration in (28): 



;^,nHn+ir 



Qni^o) 



Since n = c V^f — 1, (29) can be written as 



(29) 






2, (2/1+1) Pni^o) 



(30) 



(25) 



Figure 4 shows a plot of (30) as a function of the ratio 
(a/b) of the prolate spheroid for various (ri/6) ratios. 

The vahdity of the theoretical expression for \G 
has been corroborated by an experiment. For this 
experiment a model of the toroidal-prolate spheroid 
antenna about 2y2 in. long was constructed. The 
model was immersed in a conducting medium of salt 
water (4.89 mhos/m) and measurements of the im- 
pedance were made at the input of the transmission 
line feeding the model. Two sets of data were taken, 
the first set was over a selected frequency range for 
the toroid alone; the second set was for the toroid- 
prolate spheroid combination for the same frequency 
range. 

The electrical length effect of the transmission line 
was removed from both sets of data, then both sets of 
impedance data were converted to admittance form. 



The total tangential electric field must vanish at the 
surface of the prolate spheroid. Using this boundary 
condition, the unknown coefficients in (24) are deter- 
mined. The results are 



Cn- 



P,,(go) 



1,2, 



Substitution of Hi into (20) yields 
Fi„ = Fo + AG, 
where 



'~~'^L 



HPdl 



(26) 



(27) 




Figure 4. Graph of normalized conductance. 
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Figure 5. Graph of AG determined theoretically and experimen- 
tally. 

(a/6 = 2.06) 



A change in conductance resulting from the addition 
of the prolate spheroid to the toroid antenna was com- 
puted and is shown plotted with the theoretical AG 
in figure 5. (The effect on AG in fig. 5 caused by a 
number of turns of wire wound on the toroid was taken 
into account using the method described in the fol- 
lowing section.) 



4. Effective Length and Equivalent Circuit 

Combining the results of Swain [1965], who dis- 
cusses a toroid in a dissipative medium, and the results 
of this paper, an equivalent circuit valid for the small 
toroid-driven prolate spheroidal antenna is shown in 
figure 6. The meaning of the various circuit param- 
eters is described below: 

Ri is the self-resistance of the A'^-turn primary 
winding on the toroid. 

Li is the self-inductance of the A^-turn primary 
winding on the toroid. 

Lz is the self-inductance of the one turn secondary 
of the transformer.^ 

M is the mutual inductance between the primary 
and secondary. 
lIRm is the conductance due to the toroid by itself 
(derived by Swain [1965]). 
AG is the increase in conductance caused by the 
presence of the spheroid given by (29). 
IKojLp) is a susceptive component due to magnetic 
energy storage in the dissipative medium. 
The appendix contains a derivation of an 
approximate expression for this inductance, 
see (39). As the transformer in figure 6 
approaches an ideal transformer, the effect 
of this susceptance becomes appreciable. 



^ An investigation of the secondary resistance due to loss in the spheroid itself is not pre- 
sented here but has been performed. It showed that this resistance is negligible for typical 
practical material, e.g., spheroids made of copper and a dissipative medium of sea water or 
earth. 




Figure 6. Equivalent circuit for a toroid-driven prolate spheroidal 
antenna. 



The effective length and the driving point imped- 
ance of the toroid-driven prolate spheroidal antenna 
can be computed using elementary circuit theory once 
values have been assigned to the parameters in figure 
6. The driving-point impedance is Vi/h when / = 0, 
and the effective length is Vi/E^ when /i = 0. 

It is of interest to note that when the transformer in 
figure 6 becomes ideal, then the driving-point im- 
pedance is 



^in — ^1 + 



A^ 



and the effective length is 



4ff ~ 



Nirb^ag 



(31) 



(32) 



k^here 



Y, = llRm^AG-jl{coLj,). 



The quasi-static approximations used in this paper 
do not contain some effects that could affect signifi- 
cantly the input impedance and effective length of the 
antenna. These effects can be introduced into the 
equivalent circuit in figure 6 by the addition of appro- 
priate circuit parameters. Swain [1965] shows how 
to introduce a capacitor to correspond to parasitic 
capacities, and this paper shows how to introduce an 
inductance (Lp) to correspond to magnetic energy 
storage in the dissipative medium. Several other 
effects could be introduced depending upon the user's 
application and desire for refinement. 

5. Conclusions 

Equations (31) and (32) indicate that if the effects of 
Lp are small, and if the transformer in figure 6 ap- 
proaches the ideal, then the prolate spheroidal anten- 
na's driving-point impedance and effective length 
should be relatively insensitive with respect to fre- 
quency changes (provided the original quasi-static 
approximations are observed). A more detailed com- 
putation shows this to be so. For this computation, 
frequencies were assumed to be in the VLF band, the 
prolate spheroid was assumed to be 1 m long with an 
axes ratio (a/6) of 2.5, and the dissipative medium was 
assumed to be sea water (cr = 4 mhos/m). The sphe- 
roid was assumed to be excited by a toroid with 26 
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turns evenly distributed about it. Figures 7 and 8 
show the driving-point resistance and reactance re- 
suhing from the computations plotted versus frequency 
for several different assumed toroid permeabilities. 




12 16 20 
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Figure 7. Input resistance for computed example. 
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Figure 8. Input reactance for computed example. 



Figure 9 shows the magnitude of the effective length 
plotted versus frequency for the assumed toroid 
permeabilities. It is to be noted that a relative 
permeability of 10^ corresponds very nearly to the 
ideal transformer case. However, practical perme- 
abihties of the order of 2 X 10^ or 4 X 10^ exhibit de- 
sirable characteristics if broadband capabilities are 
desired at VLF. This example also shows that the 
Q of the antenna is 0.5 or less over the frequency 
band. 

6. Appendix 

The magnetic energy stored in the dissipative me- 
dium external to the prolate spheroid antenna is 
computed using 



W=fJH^dv, 



(33) 



where v indicates the region external to the spheroid. 
At a given frequency o), this energy can be associated 
with an inductance: 



j/2 



4cjW 



(34) 



where V is the peak voltage across the inductance. 

In principle, it is possible to substitute (25) into 
(33) and (34), and compute an expression for Lp. 
However, the computation of (33) in prolate spheroidal 
coordinates presents some numerical difficulties be- 
cause the orthogonality relations associated with 
Legendre functions cannot be used except for the 
degenerate case when the spheroid is a sphere. 

It is of some value, however, to derive formulas for 
W and Lp for the sphere, because the result can be 
extended, by approximation, to spheroids that do not 
differ too much from spheres. 

Thus, in the limit, as a prolate spheroid approaches 
a sphere, 

c->0. 




4 8 12 16 20 24 
Frequency in kilocycles 

Figure 9. Effective lengths for computed example. 



7]^ COSO, 



where c is the semifocal length of the prolate spheroid, 
and where (r, 6, cp) represents a standard spherical 
coordinate system. The result of this limiting process 
yields 



_j 

28 



Ht = HP-\-H\ 



(35) 



where 



W 



2 ^^nin + lf"^''''^'^' 
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Arr- 



2 ^^ n{n-^l) 



n^l /b\^ 



>h 



Independent computations involving prolate spheroidal 
functions performed on a digital computer show that 
the error of (39) is not too great if ajb < 3, see figure 10. 
Substituting (38) into (34), and using the fact that 
K = -V, yields 



The total magnetic energy stored external to the 
sphere is then 

W = ^\^ dr\" rdO f '" r sin d^pHl (36) 

^ Jb Jo Jo 

If the magnetic ring source is very close to the surface 
of the spheroid then n -^ b, and (36) is evaluated to be 



y, 7T(Xoa^K'b^ ^ [P/,(0)]^[l+^„]^ 



n{n+l){2n+l){2n-l) 



(37) 



This series is dominated by its first term, and the factor 
(1-h.^i) is related to the gain g given by (19). The 
series is therefore approximated to 



TTfloa^b^ 

12 



(38) 



It is now intended that (38) applies to prolate sphe- 
roidal antennas that do not depart too much from a 
sphere, and that are excited by a toroidal winding 



located 



very close to the surface of the spheroid. 




TTco^lJioO-^g^b^ 



(39) 
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Figure 10. Comparison of magnetic energy computed by two dif- 
ferent methods. 
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